Caged clusters in Aliilr4: structural transition and insulating phase 
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Aliilr4 is a phase in one of the complex metallic alloy 
systems that resist ordinary approaches to determining 
crystallographic structure or phase diagrams, because of 
two impediments to equilibration at low temperatures: 
(1) there are numerous competing complex phases of sim- 
ilar composition (2) the intermediate-temperature struc- 
tures have inherent entropy associated with block rear- 
rangements of tiles or clusters. 

Despite the measurement of 457 independent reflec- 
tions, the experimentally refined structure [T] of Aliilr4 
has twice as many sites listed as atoms, many with a fit- 
ted occupancy of less than 0.5. Until occupancy correla- 
tions are known, such structure knowledge is inadequate 
for computing electronic properties or total energies. 

Instead, in this letter wc predict the structure, aided 
by its close relationship to Al-transition metal (Al-TM) 
quasicrystals. These are described as networks of iden- 
tically oriented icosahedral clusters linked along certain 
symmetry directions [2 , e.g. the "Mackay icosahedron" 
[2] consisting of an empty center plus two concentric 
shells of full icosahedral symmetry, AI12 + Al3oMni2. 
It is generally accepted that the stable quasicrystals i- 
AlCuFe and j-AlPdMn contain a modified cluster called 
the "pseudo Mackay icosahedron" (pMI) in which the in- 
ner AI12 icosahedron is replaced by a cluster with fewer 
atoms and reduced symmetry |3j, the details of which 
are still unclear. We note that the Ca-Cd quasicrystals, 
another major class alongside the Al-TM quasicrystals, 
and complex alloys related to them, also contain icosa- 
hedral clusters with an inner shell of lower symmetry [1] 
- a tetrahedron, in that case. Thus, such inner rotatable 
clusters may be a common feature of well-ordered qua- 
sicrystal systems, and they present tractable examples 
of the above-mentioned block-rearrangement degrees of 
freedom. 

In the rest of this paper, we uncover the correct low-T 
structure of Aliilr4 using energy minimization, in several 
stages. We used both ab-initio codes (VASP package [S]) 
and empirically fitted pair potentials. First we deter- 
mined the possible atom contents of one cluster to be 
Alglr or Aliolr. Next, for each cluster, we find its atom 
arrangement and preferred orientation within the cage of 



neighboring atoms. The partially occupied sites in the 
refined structure [Tj are explained with an equal mixture 
of the two kinds of cluster. We also find the collective 
ordering of orientations and identify the ordering tran- 
sition. The equal-mixture phase is metallic, as seen in 
experiment 25J, but the composition with Alioir clusters 
is predicted to be insulating. 

As a prerequisite to exploring alternative structures 
and running extensive molecular dynamics (MD) simula- 
tions, we generated "empirical oscillating pair potentials" 
(EOFF) valid for Al-Ir in this composition range, using 
the method of Ref. [6\ to fit a database of ab-initio forces, 
calculated with VASP [5j . The database contained snap- 
shots from ab-initio molecular dynamics simulations at 
high temperatures, as well as T = structures of var- 
ious stable and unstable Al-Ir structures, in particular 
the Al2iFd8 and B2 structures. (The functional form 
and fitted coefficients of the potentials are provided in 
the supplementary material [lO'.) 

Structural optimization of individual cluster in Alulr^ 
— The reported crystal structure of Aliilr4 is a simple 
cubic arrangement of pMI clusters, each having an in- 
ner shell of ^ 40 partially occupied sites with combined 
occupancy 10. Surrounding this is the second shell, an 
Iri2 icosahedron, plus 30 Al sites slightly outside that 
icosahedron's mid-edges, of which the six sites along cu- 
bic 2-fold axes are only partially occupied. Between 
the pMI clusters, at each body center, is an Ir atom, sur- 
rounded by an icosahedron of 12 Al atoms, which are 
shared with the pMI in the role of second-shell Al. The 
key fact to determine by energy optimization is the atom 
configuration in the uncertain, inner-shell sites. 

The first question about the inner Al cluster is how 
many atoms it (optimally) contains. We addressed 
this directly by constructing initial configurations of the 
Aliilr4 unit cell with a chosen number m of inner Al 
atoms (assuming the periodicity of only one unit cell). 
The configuration was optimized, first by annealing with 
molecular dynamics (MD) using the fitted EOFF poten- 
tials at fixed volume; this was followed by ab-initio relax- 
ation using VASF [5], optimizing all structural parame- 
ters. [3. The first stage was annealing in one unit cell 
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FIG. 1: Inner clusters in Aliilr4 structures, shown as they 
order in the low-T "9.5 phase". The surrounding icosahedron 
of black atoms is the Iri2 "cage" ; another Ir atom is at center. 
The "apex" Al atoms are shown gray. (a). Alglr clusters in 
z = plane of the "9.5 phase" with their 4- fold axes oriented 
up and down along viewing z axis, (b). Alioir clusters in z = 
0.5 plane, with their 3-fold axis almost parallel to the viewing 
z axis . (c). Orientation relation of two Alglr clusters (left) 
with two Alioir clusters (right), viewed along x direction. The 
layers of identical-content clusters are stacked normal to the 
page in (a) and (b), or left to right in (c), in a checkerboard 
pattern. In the 10 phase, 10-cluster arrangement in (b) is 
repeated in 2x1x1 supercell. 



inner cluster Al 


8 


9 


9.5 


10 


11 


12 


1x1x1 


+45 


+33 




+10 


+77 


+ 132 


supercell 




+18 


+9 


+4 






flcub [A] 


7.64 


7.68 


7.70 


7.73 


7.78 


7.82 



TABLE I: Relaxed energies (in meV/atom) as a function of 
Al per pMI inner cluster. Supercell used was 2x2x2, except 
4x4x4 for the 10 Al case. Last row is mean lattice parameter, 
per fundamental cell. 



(see Table[T]), to identify gross energy differences: this ex- 
cluded the m = 8, m = 11 and to = 12 variants, leaving 
only the m = 9 and m = 10 ones which we "9-cluster" 
and "10-cluster" from here on. In a second stage, we 
annealed in a supercell to accomodate possible alterna- 
tions of the orientations, and found two nearly stable 
structures: the m — 10 filling ("10-phase") and an equal 
mixture of m = 9 and m = 10 fillings ("9.5-phase"). 
From here on, in studies of the whole structure, we limit 
ourselves to these two fillings. 

The 9-cluster is almost always a square antiprisni, with 
one square enlarged and capped by an additional Al atom 
[see FigjTJa)]. In its optimum orientation (of multi- 
plicity six), the fourfold axis is aligned with one of the 
cubic (100) axes, with the outer square of atoms in mir- 
ror planes. The 10-cluster is practically always a sort 
of trigonal prism, with one triangular end face enlarged 
and capped by an atom, and also capped on the three 
trapezoidal side faces. 



In its optimum orientation (of multiplicity 12), the 10- 
cluster is tilted rigidly so as to bring its threefold axis 
closer to (but not quite parallel to) the nearby (100) axis, 
as seen in Figure [ijb). That brings the apical capping 
Al and all three of the side-capping Al atoms (which are 
farther out from the central Ir) close to cubic (100) axes, 
and they merge with the Al atoms of the pMFs second 
shell (cage) that sit in the same direction. The pMI clus- 
ters in Fig. [l] overlap along the [100] direction such that 
some outer shell atoms of one pMI also belong to the 
inner cluster of the neighboring pMI [T^ . 

At T ^ 600K, just lOOK above the orientational Tc 
(see below), an inner Aliplr cluster spends 95% of the 
time close to one of the twelve ideal orientations de- 
scribed above (as shown by quenching a crystal using 
Alio filling in each cage.) At higher temperatures, the 
threefold axis of the 10-cluster may instead align along 
a cubic (111) axis, or a pseudo 2-fold axis (of the Iri2 
cage) . (We have not examined the high-T orientations of 
the 9-cluster.) 

The Al partial occupancies refined in [1] are in close 
agreement with our model, provided that close to half 
the clusters are 9-clusters and half are 10-clusters, and 
each is oriented independently in one of the (respectively) 
six or 12 optimum orientations. [TT] (Fig. [2]ja) shows the 
time-averaged atom density from a high— T simulation, 
containing lines of closely spaced alternative sites. The 
refined Al(3) and Al(4) sites of [T] represent this spread- 
out density, coming from four or five model sites of both 
9- and 10-clusters; the reported Al(2) and Al(5) sites 
comes respectively from one site of the 9-cluster and two 
sites of the 10-cluster in our model, showing that both 
clusters must be present. The r.m.s. discrepancy be- 
tween predicted and reported occupancies of each site is 
~12% of the occupancy (10) . 

Orientational ordering — To address the statistical 
mechanics of systems with many interacting clusters, we 
set up a replica-exchange [7] simulation in a 4 x 4 x 4 
supercell, i.e. 1920 atoms (10-phase), annealing 16 inde- 
pendent samples, each at a different temperature, span- 
ning 310K to 460K; additional single-T simulations were 
done at higher temperatures. 

The simulated alloy remains solid beyond 1700 K in 
2x2x2 supercell (the experimental melting point is above 
1600 C for xir - 0.27 [I^). The simulation lasted 4000 
cycles, each consisting of 1000 molecular dynamics (MD) 
steps with time increment At=5 is, for a total simulation 
time 20 ns. The pair-potential interaction radius was 
cut off at rc=10A. [8j. The initial state for the MD was 
our best single-cell model, repeated in all 4"^ cells. 

The low-temperature optimal structures of the 9.5- 
phase and 10-phase were determined by rapidly quench- 
ing configurations from the lowest replica (T = 310K) 
to a T = energy minimum, and selecting those with 
the lowest relaxed energies. In 9.5-phase, the structure 
was a stack along (say) the z direction of layers which 



3 



J 

\. > 

EB <» fic 

A ^ , 

> X > 





470 K 




- 480 K 




«495 K 




- 500 K 




510 K 




(b) 







-0.3 -0.296 -0.292 -0.28S 

E [eV/atom] 



FIG. 2: Disordered and ordered states of inner clusters. Phase 
transition in the "10" phase from rephca-exchange molecular 
dynamics simulation in 4 x 4 x 4 supercell of the basic (30- 
atom) unit cell, simulation results from 4x4x4 supercell of 
the basic (30-atom) unit cell, (a) time-averaged occupancy 
distributions for a slice of thickness GA taken through the 
4x4x4 supercell of the "10" phase composition, below the 
ordering transition at T =470 K (left) and above it at 510 
K (right). (b) Probability distribution P{E) for total 

energy (normalized per atom) of the "10" composition, at 
five temperatures from 470K to 510K; the double peak at 480 
K < T < 500K is diagnostic of coexistence and a first-order 
phase transition. 



structure Pearson Space parameters (A) Al Ir 

symbol group a b c 
Ref.IH cP60 P23 7.67 (cubic) 22 8 
high-T 9.5 cP178 Pm3 7.73 (cubic) 21.5 8 
9.5-phase oA236 Abm2 15.48 15.41 15.32 172 64 
10-phase oP60 Pma2 15.49 7.74 7.70 44 16 

TABLE II: Crystal data for Aliilr4 structures, from diffrac- 
tion refinement [1], the disordered high-T 9.5-state and the or- 
dered states of the 9.5-phase and 10-phase. Columns include 
lattice constants and Al or Ir atoms per cell. For "high-T" 
9.5, data in "(...)" are from Ref. [T] 



had a y/2 x a/2 checkerboard pattern of inner clusters, 
with their main (fourfold or threefold) axes alternating 
between +z and —z directions; the layers are stacked 
such that the overlaying clusters have the same orienta- 
tion of main axes. (To fully specify their orientations, 
note the 10-clusters oi +z and —z alignment are always 
related by a mirror reflection in the xy plane.) The re- 
lation between adjoining clusters is shown in Figure [l] 
In the 10-phase, the 10-clusters pair shown in Fig. [ijb) 
is repeated in the supercell doubled along y direction. 
The crystal structures are summarized in Table |TTj coor- 
dinates are available in Ref. ^TU\ . 

The ordered arrangement of the 9.5 phase contains 
alternating (square) layers of 9-clusters and 10-clusters. 
Each layer is made by repeating Fig. [ija) or (b) in both 
directions: thus, within either kind of layer, the apex di- 
rections alternate up and down. These layers are stacked 
without shifts so that chains of identically aligned clus- 
ters run perpendicular to the plane of the layers. 



We pinpointed the ordering temperature of the "10" 
phase in three ways. Firstly, we can examine the 
atom density distribution, averaged over moderate times 
(Fig. [2]Ja); for T > this has the full symmetry of the 
unit cell, but for T < this shows (despite some fluc- 
tuations) a clear symmetry breaking (to the cell-doubled 
structure just described). A second evidence for a phase 
transition is available within the replica-exchange simu- 
lation: a larger energy spacing A-B; = E{Ti) — -B(Tj_i) 
between the replicas at consecutive temperatures, sug- 
gesting a latent heat. The third and most convincing 
signature is that, for a flnite system, the ensemble at 
temperatures close to Tc is a mixture of the two phases, 
weighted according to the difference in their free ener- 
gies. This is evident in Figj2|b), from which we can read 
off Tc < 495-fsr. The peak separation in Figj2jb) shows 
that the ordered and disordered states in the 10-phase 
differ by (a latent heat of) 3 meV per atom, requiring an 
entropy difference at Tc of (30) (3 meV)/Tc sa ln(9) per 
30-atom cell. That is, we have effectively Ncs ~ 9 states 
per cluster, comparable with Ncs = 12 ideal orientations. 

The 9.5 phase has two kinds of order - the 9A1/10A1 al- 
ternation in cluster content and the orientations - which 
might appear in separate transitions. The first order- 
ing requires vastly longer equilibration times (for inter- 
cluster Al diffusion) and we were unable to identify any 
sharp transitions. 

Phase stability in Al-Ir system — We now turn 
to the T = phase stability of Aliilr4 and other Al- 
Ir compounds. The currently accepted Al-Ir phase di- 
agram |13j with < xir < 0.5 shows six compounds as 
stable: Al9lr2 (in the AlgCo2 mP22 structure); Allr; and 
around xit- w 1/4, there is Alair plus the complex phases 
Aliilr4, Al28lr9 [14 and orthorhombic Al45lri3 [TS]. We 
additionally tried other complex phases structures not re- 
ported in Al-Ir, among them : (i) Al2ilr8 in the Al2iPd8 
structure, [16] which is a packing of sixteen of our 10- 
clusters per unit cell, without any icosahedral cages (and 
sharing a few Al atoms), (ii) Al4ilr23 in the Al4iCu8lri5 
structure [TH], equivalent to Al68Pd2oRui2 [H] and Al- 
CuRuSi [20]; here the even- vertex clusters alternate be- 
tween an Alioir trigonal cluster and an AlCu8Al6 cluster, 
identical to atomic arrangement in B2 cubic structure. 

Of these, the available structures of Aliilr4 and Al28lr9 
include many fractional sites, so in total energy calcu- 
lations we had to try versions of these structures with 
various ways of realizing the occupations. The relaxed 
total energies of all these compounds (and the pure el- 
ements) were computed using VASP [S]. We found 
Allr and Alglr2 to be stable, as expected; the only stable 
phases around xir ~ 1/4 were A^slria and (surprisingly) 
Al2ilr8. Relative to a corrected tie-line including Al2ilr8, 
we found Alsir to be unstable by 25 meV/atom; this is re- 
duced to only '-^ 8 meV/atom in a variant (Ali7lrg) with 
a tripled cell in which one Al is removed. Thus, in contra- 
diction to Ref. QBJ we believe Alsir is a high-temperature 
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FIG. 3: Electronic density of states (DOS) of T = or- 
dered structures, where _E = is the Fermi energy; (curves 
are smoothed using the Methfessel-Paxton smearing method). 
The 10-phase has a gap of ~ 40 meV (blown up in inset, tetra- 
hedron method) but in the 9.5 phase, the Fermi level is ~ 0.5 
eV shifted from a pseudogap. 

phase only. The phases Al4ilr23, Al28lr9, and Aliilr4 in 
the 9.5-phase were unstable by rather small amounts of 
6, 7, and 9 meV/atom, respectively; in the latter two, 
the site disorder contributes entropy which may explain 
their stability at higher temperatures. The heat of for- 
mation for Aliilr4 is large: —0.738 eV/atom; for Al2ilrg 
it is —0.760 eV/atom (sec Supplementary material for 
detailed information on formation energies and stability 
data for all discussed phases). 

Electronic density of states — A specific interest in Al~ 
Ir and certain other Al-TM systems is the possibility of 
an insulating alloy, the elemental constituents of which 
are all good metals, e.g. AI2RU in the TiSi2 structure [21; 
Al2Fe in the MoSi2 type-structure ^21,,22J. Quasicrystal 
i-AlPdRe (built from pMI clusters) is long claimed to be 
a semiconductor p3l I24j . Predicting gap formation in 
these alloys depends critically on the accurate relaxation 
of atomic positions |23) . 

We find the Aliilr4 electronic density of states (DOS) 
in the ordered (T = 0) "10" phase has a narrow gap, 
and significant pseudogap as shown in Fig. [3j so this 
alloy should be a semiconductor. By contrast, the or- 
dered "9.5" phase is predicted to be robustly metal- 
lic [Fig. [sjjb)] (In the disordered high-temperature "10" 
phase [not shown] the gap tends to get filled in, so that 
phase is also likely to be metallic.) We call attention to 
the fact that these very similar structures have radically 
different electronic properties, all due to details of the 
placement of certain Al atoms that link adjacent clus- 
ters. 

Discussion — In conclusion, starting from diffraction- 
data refined average structure of AliiIr4[T], we com- 



bined molecular dynamics simulations with pair poten- 
tials to discover a well-defined, asymmetric inner Alglr 
and Alioir clusters, with variable orientations, and an 
ordering transition to new orthorhombic phases. We sug- 
gest an experimental search for the ordered phase having 
all clusters of Alioir type, as it has a gap exactly at the 
Fermi energy, and hence should be a semiconductor. We 
also found that, contrary to the accepted phase diagram, 
the stable low-temperature phase at xij. 1/4 is Al2ilr8. 

Several other known structures are made by placing 
pMI clusters and/or variants on the same simple cu- 
bic lattice with a ~7.7A: not only AUiCuglris, [18], 
Al68Pd2oRui2 [12 and AlCuRuSi [2D], which we men- 
tioned, but also Al7oPdioFe2o [55], and e-Ag7Mg26 (in 
the latter case, the "even" -node inner clusters are AgMgg 
cube, while "odd" nodes are perfect Mackay icosahe- 
dra). 

Furthermore, the lowest energy version yet found of 
the i-AlPdMn quasicrystal structure [5BJ [37] consists of 
exactly the same pMI clusters described above, in partic- 
ular (i) they have a short (~ SA) linkage along the icosa- 
hedral twofold direction, rather than the ~ 12. sA linkage 
known from a-AlMnSi [2] (ii) they have an AligMn inner 
cluster in the same shape as the one described here, which 
appear in various orientations with only small energy dif- 
ferences. The decagonal-related ^'(AlPdMn) phase con- 
tains pMI clusters with AlgMn and AlipMn inner clus- 
ters, not quite identical to those we presented here. 

The icosahedral cage containing an inner cluster that 
breaks icosahedral symmetry is known in a quite different 
family of quasicrystals and related compounds: CaCdg 
or ScZne, with an inner Cd4 or Zn4 tetrahedron [^ I30|. 
Those, too, show orientational orderings at low temper- 
atures [3TI [32]. We observe that, if there is any way at 
all in icosahedral quasicrystals to implement the long- 
sought "local matching rules" that would stabilize an 
ideal quasiperiodic ground state, or merely rules that 
limit tiling randomness to large spatial scales the likeliest 
candidate are these inner clusters that spoil the clusters' 
high symmetry. 

On the theoretical side, this work on Al-Ir could be ex- 
tended in the following directions: (1) simulation of the 
inner-cluster dynamics for comparison to neutron data, 
analogous to the Cd4 tetrahedron in CaCde [3^; (2) fit- 
ting an orientational interaction between adjacent clus- 
ters [33] for Monte Carlo estimation of the orientational 
Tc- (3) a possible Al-Ir quasicrystal based on these clus- 
ters. Experimental clarification of the low-temperature 
stabilities of Al-Ir around xjr ^ 0.27 are of course highly 
desirable. 
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tural data on Al4iCu8lri5 phase prior to publication, R. 
Tamura for fruitful discussions, and Yu. Grin for care- 
ful reading of the manuscript and comments, in partic- 
ular pointing out cluster relationship with AlgMn struc- 
ture. This work was supported by DOE Grant DE-FG02- 
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SUPPLEMENTARY INFORMATION 

In this section, we give tables of three kinds of detailed 
information omitted from the main text of the paper. 



Al-Ir phases for stability calculations 

Table [III] gives the data used to find the predicted bi- 
nary phase diagram of the Al~Ir. system, by the usual 
convex-hull construction. Despite the structure names, a 
binary Al-Ir composition is used in all cases, with Ir re- 
placing any transition metals. All energies are ab-initio 
calculations with the VASP code. 



structure 


AE 


AH 


x(Al) 


x(Ir) 




meV/at. 


meV/at. 


% 


% 


Al9Co2.mP22 





-544.0 


81.8 


18.2 


Al45lri3.oP236 





-658.8 


77.6 


22.4 


Al2iPd8.tI116 





-759.6 


72.4 


27.6 


AlIr.cP2 





-955.3 


50.0 


50.0 


Al5Co2.hP28 


2.8 


-765.4 


71.4 


28.6 


AliiIr4.oP60'' 


3.7 


-737.9 


73.3 


26.7 


Al4iIr23.hR64 


5.6 


-826.9 


64.1 


35.9 


Al28lr9.hP236 


7.2 


-672.3 


76.5 


23.5 


Al3lr.hP24'' 


7.8 


-722.5 


73.9 


26.1 


AliiIr4.oA236'= 


9.2 


-741.2 


72.9 


27.1 


Ablr.hPS 


25.0 


-684.2 


75.0 


25.0 



" "10-phase" 

triple supercell of Al3lr, with one less Al atom (site vacant). 
^ "9.5-phase" 

TABLE IH: Energy diagram of Al-Ir system at T=0K. Col- 
umn AE is energy in meV/atom by which a structure is un- 
stable relative to a mixture of competing stable compounds. 
Stable compounds forming convex hull of energy-composition 
scatter plot have AE — 0. 
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Ci 




Vi 


C2 ri2 


fc* 


Al- 


^Al 


9.9811X 


10^ 


8.7044 


-11.3044 5.4518 


3.5646 2.5119 


AI 


-Ir 


1.3764X 


10^ 


13.3271 


9.2649 3.8320 


3.1919 1.1091 


Ir- 


-Ir 


3.2083X 


10^ 


12.7731 


-7.2818 3.0765 


-2.8274 0.4493 



TABLE IV: Fitted parameters for Al-Ir EOPP potentials. 



Fitted pair potentials 

Most calculations in the paper were based on em "em- 
pirical oscillating pair potentials" (EOPP) of form 



Vir) 



C2 



cos(fc*r -t- 0* 



(1) 



as presented by Mihalkovic and Henley (main text Ref. 
[6]). where r is the distance between a pair of atoms. 



The fitted coefficients we used are listed in Table IV Fig. 
|4] gives a scatter plot demonstrating the goodness of fit 
and a plot of the potentials themselves. As is typical of 
Al-transition metal (Al-TM) potentials, with TM=Ir in 
this case: (a) the Al-Al potential has no nearest-neighbor 
well but only a shoulder (b) the Al-Ir potential has a 
very deep well at the nearest-neighbor distance (c) the 
Ir-Ir potential is unfavorable for nearest neighbors but 
has a deep well at the second neighbor distance (d) all 
the potentials have relatively strong Friedel oscillations. 



Crystal structure of phase with Alioir clusters 

Table [V] lists all the Wyckoff positions of the low- 
temperature structure of the "10-phase" of Aliilr4; The 
cluster orientations correspond to Figure 1(b) of the main 



0.2 





— Al-Al (pot) 

— Ir-Al(pot) 

— Ir-Ir (pot) 

A 



site 
Irl/2b 
Ir3/2c 
Ir5/4d 
All/4d 
A13/4d 
A15/4d 
A17/4d 
A19/4d 
Alll/4d 
A113/2C 



X 
0.5000 
0.2500 
0.5028 
0.5006 
0.3521 
0.5898 
0.3668 
0.4059 
0.6644 
0.2500 



0.5000 
0.4745 
0.7269 
0.7884 
0.6752 
0.4900 
0.8794 
0.9910 
0.7314 
0.6492 



0.5135 
0.2897 
0.0113 
0.6822 
0.5022 
0.8229 
0.2035 
0.8849 
0.0690 
0.0134 



site 
Ir2/2c 
Ir4/2c 
Ir6/4d 
A12/4d 
AI4/4d 
AI6/4d 
A18/2C 
A110/2C 
A112/2C 



X 
0.2500 
0.2500 
0.3868 
0.3544 
0.5850 
0.5175 
0.2500 
0.7500 
0.2500 



0.9952 
0.5099 
0.9996 
0.3333 
0.4755 
0.8086 
0.1771 
0.1536 
0.0944 



0.0140 
0.7156 
0.5274 
0.5161 
0.2149 
0.3328 
0.7409 
0.7156 
0.3570 



TABLE V: List of Wyckoff sites for "10-phase" structure. 
Space group is Pma2 (#28), Pearson symbol 0P6O. VASP- 
optimized lattice parameters are a= I5.492A, 6=7.725Aand 
c=7.685A. 



site 


X 




Y 


Z 


site 


X 




Y 


Z 


Irl/8d 


0.7498 





0006 


0.0008 


Ir2/4c 


0.0003 





2500 


0.7498 


Ir3/4c 


0.7499 





2500 


0.6054 


Ir4/8d 


0.3602 





0013 


0.7517 


Ir5/8d 


0.5001 





8857 


0.4935 


Ir6/4c 


0.7506 





7500 


0.4015 


Ir7/4c 


0.5018 





7500 


0.2503 


Ir8/4c 


0.2619 





2500 


0.6151 


Ir9/4c 


0.2461 





7500 


0.3962 


Irl0/8d 


0.8581 





9789 


0.2506 


Irll/8d 


0.0006 





3900 


0.4985 


All/8d 


0.3968 





0026 


0.9159 


A12/8d 


0.1609 





3489 


0.5097 


A13/8d 


0.2589 





0888 


0.1491 


A14/8d 


0.2355 





9163 


0.1534 


A15/8d 


0.1613 





6470 


0.4972 


A16/8d 


0.2433 





0876 


0.8458 


A17/8d 


0.3353 





6469 


0.5003 


A18/8d 


0.3340 





3545 


0.5049 


A19/8d 


0.2591 





9076 


0.8499 


A110/8d 


0.4067 





9853 


0.0921 


Alll/8d 


0.9012 





8370 


0.3482 


A112/4C 


0.5471 





7500 


0.5829 


A113/8d 


0.4266 





8563 


0.6523 


A114/4C 


0.4199 





2500 


0.3969 


A115/4C 


0.5860 





2500 


0.3900 


A116/8d 


0.4944 





9075 


0.8078 


A117/4C 


0.3149 





2500 


0.2534 


A118/8d 


0.9035 





9986 


0.9101 


A119/8d 


0.9061 





9990 


0.0891 


A120/8d 


0.3622 





8352 


0.2277 


A121/8d 


0.0296 





9100 


0.2485 


A122/4C 


0.0351 





2500 


0.5873 


A123/8d 


0.8927 





1643 


0.6627 


A124/4C 


0.1662 





7500 


0.2510 


A125/4C 


0.0534 





7500 


0.4062 



TABLE VI: Atomic structure of "9.5-phase" (oA236) in 
space group Abm2 (# 39). Optimized lattice parameters are 
a=15.321 A, 6=15.475 A, c=15.412 A. Cluster centers are Ir2 
( "9-clusters" ) and Ir7 ("lO-clusters"). 
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FIG. 4: EOPP potentials fitted to ab-initio (VASP) force 
and energy datapoints. Fit to forces (left panel) is shown 
with fitted forces (vertical axis) plotted against VASP forces 
(horizontal axis). The fitted pair potentials are shown in the 
right panel. Parameters 



text, corresponding to a cubic cell doubled in one direc- 
tion (becoming the a direction) ; since the cubic structure 
lacks 4-fold symmetry, this causes the b and c cell param- 
eters to become different by ~ 0.6%. 



site X Y Z 
Irl/4b 0.5000 0.0000 0.2500 
Ir3/8j 0.9954 0.3955 0.5000 
Ir5/16k 0.2503 0.5011 0.3603 
Ir7/4a 0.0000 0.0000 0.7500 
A12/16k 0.0877 0.9954 0.3987 
A14/16k 0.9124 0.5012 0.1026 
A16/16k 0.1544 0.0893 0.7501 
A18/16k 0.1403 0.3379 0.4147 
A110/8j 0.2925 0.4107 0.5000 
A112/8j 0.4193 0.2521 0.5000 



site X Y Z 
Ir2/8j 0.2497 0.2507 0.5000 
Ir4/16k 0.3926 0.2502 0.7704 
Ir6/8j 0.9997 0.1036 0.5000 
All/16k 0.8442 0.5900 0.2517 
A13/16k 0.4831 0.8460 0.1605 
A15/16k 0.9870 0.1541 0.3407 
A17/16k 0.1490 0.1544 0.4135 
A19/16k 0.2780 0.2536 0.6603 
Alll/8j 0.3035 0.0941 0.5000 



TABLE VII: "9-phase" structure (o/232), space group Ibam 
(# 72), orthorhombic lattice parameters are o=15.136 A, 
6=15.517 A, c=15.457 A. "9-cluster" centers are at Ir2 site. 



site X Y Z 
Irl/3a 0.0000 0.0000 0.0020 
Ir3/3a 0.0000 0.0000 0.2496 
Ir5/9b 0.6312 0.0770 0.4618 
Ir7/9b 0.3639 0.9178 0.5390 
Ir9/9b 0.7434 0.8712 0.4665 
Irll/9b 0.2580 0.1292 0.5313 
A12/9b 0.8717 0.1251 0.4363 
A14/9b 0.2697 0.2245 0.2304 
A16/9b 0.0192 0.1585 0.1729 
A18/9b 0.7519 0.7866 0.7716 
A110/9b 0.9820 0.8361 0.8292 
A112/9b 0.5229 0.1744 0.4054 
A114/9b 0.8615 0.5862 0.3442 



site X Y Z 
Ir2/3a 0.0000 0.0000 0.7488 
Ir4/9b 0.2685 0.2430 0.6356 
Ir6/9b 0.2696 0.9981 0.3618 
Ir8/3a 0.0000 0.0000 0.5959 
Irl0/3a 0.0000 0.0000 0.4027 
All/3a 0.0000 0.0000 0.4995 
A13/9b 0.1273 0.2493 0.5613 
A15/9b 0.0399 0.2553 0.2678 
A17/9b 0.1447 0.1397 0.3285 
A19/9b 0.9571 0.7422 0.7301 
Alll/9b 0.8418 0.8617 0.6704 
A113/9b 0.1543 0.4314 0.6314 
A115/3a 0.0000 0.0000 0.9080 



TABLE VIII: Structure of rhombohedral Al4ilr23 {hR64) 
phase, space group R3 (# 146) in hexagonal setting. Op- 
timized lattice parameters arc a= 11.032 Aand c=27.073 A. 
"10-clusters" surround Irl sites. All site is center of"B2"- 
structure local environments. 



